Adiabatic cooling of trapped nonneutral plasmas 
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Nonneutral plasmas can be trapped for long times by means of combined electric and magnetic 
fields. Adiabatic cooling is achieved by slowly decreasing the trapping frequency and letting the 
plasma occupy a larger volume. We develop a fully kinetic time-dependent theory of adiabatic 
cooling for plasmas trapped in a one-dimensional well. This approach is further extended to three 
dimensions and applied to the cooling of antiproton plasmas, showing excellent agreement with 
recent experiments [G. Gabrielse et al., Phys. Rev. Lett. 106, 073002 (2011)]. 
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Introduction. — Nonneutral plasmas (consisting purely 
of electrons, ions, or positrons) have received much at- 
tention in the last few decades 1], particularly because 
very accurate table-top experiments can be devised to 
study their properties. Nonneutral plasmas can be con- 
fined for much longer periods of time (several days) com- 
pared to quasi-neutral plasmas and are more easily cooled 
down to low temperatures. Their confinement is usu- 
ally achieved by means of a cylindrical device (Penning- 
Malmberg trap) , in which an axial magnetic field ensures 
the radial confinement while the axial trapping is pro- 
vided by a set of appropriately distributed electrodes. 

The creation and confinement of nonneutral plasmas 
consisting of positrons or antiprotons (p) is a major is- 
sue for the production of significant amounts of neutral 
antimatter in the form of antihydrogen (H) atoms. Spec- 
tacular advances in this field have been achieved in recent 
years by several international collaborations, such as AL- 
PHA 0, ATRAP U, and ASACUSA [4]. The motiva- 
tion of this recent and exciting research is to compare the 
physical properties of matter and antimatter, for instance 
the spectral lines of antihydrogen or its behavior under 
the influence of gravity. The recently established inter- 
national collaboration GBAR aims at determining the 
gravitational acceleration of H atoms byletting them fall 
in the gravitational field of the Earth Q. For the same 
purpose, the AEGIS collaboration proposes to measure 
the shift of the interference pattern of free-falling antihy- 
drogen atoms going through a series of gratings @. 

The H atoms are usually trapped in a magnetic well, 
where the effective confining potential has a depth of 
about IK in temperature units Q. As the atoms are 
created with much higher energies, their cooling consti- 
tutes a major challenge to achieve an effective confine- 
ment. The ALPHA collaboration Q recently reported 
confining several hundred H atoms for more than 15 min- 
utes. The antiprotons were first pre-cooled by interacting 
with cold electrons and then brought down to 40 K by 
evaporative cooling Q. Unfortunately, during evapora- 
tive cooling most charged particles are lost, leading to 



a lower p density. In contrast, using an adiabatic cool- 
ing technique, Gabrielse et al. [8| managed to cool down 
3 x 10 6 antiprotons to about 3K, with almost no losses 
observed. 

Adiabatic cooling is achieved by lowering the trap- 
ping frequency uj that confines a nonneutral plasma 
within a one-dimensional (ID) harmonic well, U{z) — 
muj 2 (t)z 2 /2. If the plasma density is low enough, one 
can adopt a single-particle approach. For a slowly vary- 
ing frequency, the action E/ui is conserved, so that the 
energy (and the temperature) of the plasma decreases lin- 
early with the frequency, i.e., T ~ uj. When the plasma 
density is high enough, the effect of the space charge can- 
not be neglected and the single-particle approach breaks 
down. Li et al. Q have developed a thermodynamic 
theory of adiabatic cooling for an ion cloud in a Penning 
trap, assuming that the ions always remain at Maxwell- 
Boltzmann equilibrium [10| . 

The purpose of this Letter is to develop a fully kinetic 
description of the adiabatic expansion and cooling of a 
nonneutral plasma, without making any assumptions on 
its thermodynamic properties except for the initial equi- 
librium. The results will be compared to recent experi- 
ments from the ATRAP collaboration [8j . 

Kinetic model in ID. — In most experiments, the 
charged particles are confined in a harmonic potential 
well, with a strong magnetic field along the axial (par- 
allel) direction z, which ensures radial (perpendicular) 
confinement. It is then reasonable to restrict our analysis 
to a ID geometry along the z axis. Nevertheless, particle 
collisions can couple the axial and radial motion, lead- 
ing to energy equipartition between the parallel and per- 
pendicular temperatures. Thus, from a thermodynamic 
point of view, the expansion is ID when collisions are 
rare, and 3D when collisions are dominant 

We first concentrate on the collisionless ID regime, 
where the antiproton plasma can be described by the 
Vlasov-Poisson equations: 
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where — e and m are the antiproton charge and mass, 
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E z is the parallel self-consistent electric field, a±_ is the 
2D number density in the perpendicular plane, and n = 
J fdv z is the ID density along z. The initial condition is 
computed self-consistently by assuming that the plasma 
is at thermal equilibrium with temperature Tq. 

If the frequency lo decreases in time, the plasma will 
experience an expansion. In order to follow this expan- 
sion, it is useful to use the scaling techniques introduced 
inRefs. [TUG]!- We define the new phase space variables 
(£, 77) and the new time 9 in the following way: 

z = C(t)£; dt = A 2 (t)d8; v z = (C/A 2 )r] + C£. (3) 

Our aim is to find the "right" scaling transformations 
such that the plasma will be frozen in the scaled vari- 
ables. Then, the scaling factor C(t) will inform us on 
the expansion law, which in turn can be related to the 
evolution of the plasma temperature. 

We also need to transform the distribution function 
in order to guarantee the invariance of the total number 
of particles, f(z,v z ,t)dzdv z = F(£,n 1 9)d£dri. Finally, 
the scaled electric field is defined as £(£,,0) = E z (z,t) 
in order to keep Poisson's equation unchanged: d^£ — 
~ S J J " I Fdr). In the scaled variables, the Vlasov equation 
reacfs as follows: 

W +7? f" + ^ [O(e ' ?7 ' W? ' ?7 ' 0)]=0 ' W 
where the scaled acceleration is 
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(5) 

Let us now introduce some specific expressions for the 
scaling coefficients: A(t) = (1 + flt) a and C{t) = (1 + 
ilt) s . In addition, we postulate a similar form for the 
time-dependent trapping frequency: uj(t) = aj (l+^)~' 3 , 
where a, /3, and 6 are real positive numbers. Here, fi > 
is a parameter that determines the rapidity of the change 
of the trapping frequency. With these prescriptions, the 
scaled acceleration becomes: 

Q = -(l + nt) 4a - s e£/m 

- [6(6 - i)n 2 (i + nt) 4a - 2 + + nt) 4a - 2 ^ 

+ 2n(a-5)(l + nt) 2a - 1 T 1 . (6) 

The idea is to choose the parameters a and 6 so that 
as many terms as possible in Eq. ^ become time- 
independent. For all values of j3, we need to choose 
a = 5/4 so that the coefficient of the first term is con- 
stant. Then, two cases clearly appear. For j3 > 1 one 
should take 6 = 2, whereas for /3 < 1 one should have 
6 = 2/3. Note that for j3 = 1 all coefficients become 
constant. 

The transformed equations describe a "virtual" physi- 
cal system made of charged particles interacting through 
(i) the self-consistent electric field £, (ii) a harmonic field 
composed of a constant and a decreasing-in-time term, 



and (iii) a friction field, either constant or decreasing 
with time. In the scaled phase space, such a system will 
relax to a steady state with zero velocity (because of 
the friction) and uniform spatial density uq (defined by 
e 2 a±no/mso = 0Jq), thus creating a linear repulsive force 
£ that exactly cancels the second term of Eq. ^ . 

The approach to the asymptotic solution is 
nicely seen in Fig. 1, where we plot the initial 
and final densities. The asymptotic profile is dic- 
tated by the structure of the scaled equations ([4])- 
(|6j) and thus constitutes a universal attractor to- 
wards which any initial condition (not necessar- 
ily an equilibrium) should converge. This is true 
even for rapid variations of u}(t), thus generaliz- 
ing results that rely on assumptions of thermal 
equilibrium and adiabatic expansion [9|]. 
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FIG. 1: Density profiles in the scaled space at the beginning 
(dashed line) and at the end (solid line) of the simulation. 

As the system becomes frozen in the £-space, the ex- 
pansion law is determined by the scaling function C(t) 
and hence the exponent 6: z ~ t 2 @ for (3 < 1 and z ~ t 2 
for P > 1. It is easy to understand why the value (3 = 1 
plays a pivotal role: the expansion is adiabatic when 
<C UJ, which implies 0CI < u (l + Qt) 1 " 13 . When 
/3 > 1, the previous expression cannot be satisfied for all 
times. In contrast, for /3 < 1, the expansions is always 
adiabatic, provided that j3Q < loq. In the following we 
will concentrate on such adiabatic regimes. 

In the scaled space, the density approaches asymptot- 
ically the value n(£) = no Using the conservation of 
the total number of particles n{z)dz = n(£)d£ to trans- 
form back to real variables, one finds that the density 
behaves as n(z,t) = no(l + flt)~ 2 P . Finally, using the 
isentropic law for a weakly coupled plasma T ~ n 7_1 
[with 7 = (d + 2)/d = 3 in ID] and the expression for 
u>(t), one obtains the temperature as a function of the 
frequency: T/Tq = (lo/luq) 4 . 

An important consequence is that the plasma fre- 
quency behaves as: uj p (t) = (e 2 f7j_n(i)/me ) 1 ^ 2 = 
Wo/(l + Qt) 13 = ^(t). In other words, the plasma fre- 
quency approaches the harmonic frequency and remains 
locked to it for the rest of the evolution. 

The relation T ~ w 4 represents an asymptotic law. 
However, some transient regime must exist in the case of 
low-density plasmas, for which the Coulomb interaction 
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is almost negligible. In this case, the expansion should 
initially follow the single-particle law T ~ ui. In contrast, 
for high-density plasmas the asymptotic law should be 
observed from the beginning of the evolution. These two 
regimes can be characterized by the total number of par- 
ticles in ID, N — J fdzdv z (normalized to noAo, where 
A = y/k B T /m/uj ). 

In order to check the above conjectures, we performed 
numerical simulations of the nonlinear Vlasov-Poisson 
equations (P)-©, starting from self-consistent equilibria 
characterized by different values of N. The equations 
were solved in the scaled phase space (£, 77) using an Eu- 
lerian Vlasov code [l5j . 

In Fig. 2, we plot the temperature at the center of the 
plasma cloud (z = 0) against the inverse harmonic fre- 
quency, for (3 = 1, Q, = 0.02w , and three values of N. As 
expected, for N < 1 the temperature evolves first as u> 
(single-particle transient) and later as u 4 (self-consistent 
asymptotic behavior). For N > 1, the asymptotic behav- 
ior is observed right from the start. The same pattern 
was observed for runs with j3 < 1. 
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FIG. 2: Plasma temperature as a function of the inverse of 
the harmonic frequency cj in the ID regime, for N — 10 (a), 
N = 0.3 (b), and N = 0.1 (c). The dashed straight lines have 
slopes —1 and —4. 

Extension to 3D. — The above model, Eqs. (HI)-®, is a 
mean-field approximation that neglects binary collisions 
between particles. However, collisions are important in 
antiproton cooling experiments and their main effect is 
to drive the system towards energy equipartition between 
the parallel and perpendicular motion. In order to model 
collisions, we add a relaxation term on the right-hand side 
ofEq. Q: 
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where u(z,t) = J fv z dv z /n, T±(t) is the perpendicular 
temperature, which obeys the following equation: 



dT 



(S) 



and ksTAt) = mj /(0, v z ,t)v z dv z /n(0,t) is the parallel 
temperature at the center of the plasma. A factor 1/2 
appears in Eq. ©, but not in Eq. ([7]), because the 



perpendicular motion is associated with two degrees of 
freedom. 

In ID, the temperature depends on the frequency as 
T ~ w 4 and the density as 11 ~ a; 2 , implying that 
Lo p (t) ~ u>(t). This relation does not depend on the di- 
mensionality of the system and can be extended to the 
3D case. Using the isentropic law in three dimensions, 
T ~ ?i 2 / 3 , yields the expression T ~ w 4 / 3 . This is the law 
that we expect to observe in the collisional 3D regime. 

In order to check this conjecture, we solve Eqs. ([lj, 
([7]) and ||5J) using the same scaling techniques in- 
troduced above. Notice that in Eq. © we have two 
free parameters, namely a and 7. The latter determines 
the expansion of the plasma density and should be kept 
the same as in the ID simulations, i.e. 7 = 2/3 for 
the adiabatic regime. The exponent a can be chosen 
so that, in the scaled variables, the distribution func- 
tion is bounded in velocity space. The relationship be- 
tween the temperature in the real and scaled spaces is 
as follows: T(£,0) = (A 4 /C 2 )T{z, t). If T ~ w 4 / 3 , then, 
in order for the scaled temperature T to remain con- 
stant, we need to take a = 4/3/3. In addition, by choos- 
ing (3 = 3/8 several terms in Eq. © are algebraically 
simpler, for instance the friction term becomes time- 
independent. With these parameters, the distribution 
function should remain bounded in a finite phase-space 
volume in the scaled variables, thus greatly simplifying 
the numerical integration without any loss of generality 
as to the physics. Equations (J7J)-((U) are also scaled using 
the same technique. 

The numerical results confirm our conjecture. The 
case depicted in Fig. 3 represents a simulation with 
f2 = 0.02wo, N — 30, and u co a — 0.001w - Initially, when 
the collision term has not yet had time to act, the par- 
allel temperature follows the ID law T« ~ ui 4 , while T± 
remans almost constant. On a longer timescale, however, 
the collisions drive the system towards energy equiparti- 
tion, and both temperatures start decreasing as w 4 / 3 . For 
larger collision rates, the law T ~ w 4 / 3 is observed from 
the beginning of the simulation. The total temperature 



is defined as T 



T, 



3^ X . 

Comparison to experiments. — We now have all the el- 
ements to perform a detailed comparison with the recent 
experimental results of Ref. @, where Np = 5 x 10 5 an- 
tiprotons were confined in a cylindrical trap of volume 
V = nR 2 L z , with radius R w 2 mm and length L z m 10 
mm. The axial magnetic field was B = 3.7 T. Each 
point on Fig. 3 of Ref. Q corresponds to a different 
experiment in which the harmonic trap frequency is low- 
ered from ujq/2tt (which ranges bewteen 90 kHz and 3 
MHz) to the final value cj/27t = 75 kHz. During the en- 
suing expansion, the antiproton plasma cools down from 
To = 31 K to a final temperature T. The total duration 
of the cooling is tf = 100 ms. 

The experiments show that the final temperature be- 
haves as T ~ lj 1 - 2 , which is rather close to the law 
T - w 4 / 3 that arises from our 3D model. In addition, 
the final temperature levels off at a value T w 3.5 K 
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FIG. 3: (Color online). Parallel, perpendicular, and total 
temperatures as a function of the inverse of the instantaneous 
trapping frequency to(t) in the 3D regime. The straight lines 
have slopes —4 and —4/3. 

and does not decrease any further for frequencies above 
ujq j '2tt = 500 kHz. This unexpected saturation is still 
unexplained and was tentatively attributed by the au- 
thors to intrinsic limits on the measuring apparatus, to 
some source of technical noise, or to some unspecified 
new physics appearing at high frequencies. 

In order to perform our simulations, we need to spec- 
ify three dimensionless parameters, namely N, Q/ujq, 
and v co \\/uj a . In 3D, TV = f fdzdv z is related to the 
total number of antiprotons by the expression: TV = 
Np/(no\oa±irR 2 ). The parameter SI, which determines 
the rapidity with which the harmonic frequency de- 
creases, can be expressed in terms of the cooling time 
tf by using the definition of ui{t). 

Finally, the collisional equipartition rate was deter- 
mined using the detailed calculations of Ref. [l6j]. For 
an unmagnetized plasma, the rate is roughly i-Wmag = 
v / 2k B T/m(N p /V)b, where b = e 2 /(2Tre a k B T), which 
yields ^ unm ag ~ 3300 s _1 for the experimental param- 
eters of interest here. The effect of the magnetic field 
is to reduce substantially the collision frequency. For 
B = 3.7 T, the correction factor can be estimated to be 
around 0.2 [l6j], yielding v co \\ rj 660 s _1 . 

To compare with the experimental results, we per- 
formed a set of runs for different values of luq , which cor- 
respond to different dimensionless parameters. Two such 
runs are displayed as solid lines on Fig. 4. In the same 
figure, we show the final temperature after adiabatic cool- 
ing as a function of the initial trapping frequency. The 
agreement between the simulations and the experimental 
results from Ref. [§| is excellent. Not only is the slope 
of the temperature curve T(u!q) correct for small values 
of the harmonic frequency, but the leveling off at high 
frequencies is also reproduced with great accuracy. 

Physically, this saturation is linked to the strength of 



the collision rate relative to ujq. On the left of the di- 
agram, v co \\/ljo is relatively large, so that equipartition 
occurs early and the corresponding points are roughly 
aligned along a line with slope 4/3 (see the curve for 
u /2tt = 470 kHz on Fig. 4). On the right of the dia- 
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FIG. 4: (Color online). Final temperature T after adiabatic 
cooling against the initial harmonic frequency ujq/2-k. Black 
squares: present simulations; open diamonds: experimental 
results from The red triangle corresponds to a case where 
v co \\ has been increased by a factor 5. The thin solid lines 
show the temperature evolution as a function of uq/ui (top 
scale) for two values of wo. 



gram, v co \\/u>q is smaller and the temperature undergoes 
a transient with shallower slope (see curve for ojq/2tt = 2 
MHz on Fig. 4 and also the total temperature curve 
on Fig. 3), during which the parallel and perpendicular 
temperatures have not yet equilibrated. This results in 
a final total temperature that is higher than what could 
be expected from the law T ~ w 4 / 3 . Indeed, by artifi- 
cially increasing u co \\ by a factor 5, we could lower the 
final temperature to roughly 2 K (red triangle in Fig. 
4). Since the collision rate is very sensitive to the exter- 
nal magnetic field, this suggests that lower temperatures 
may be achieved by reducing the latter, a conjecture that 
could easily be tested experimentally. 

In summary, we have presented a fully kinetic theory 
for the adiabatic cooling of charged particles both in ID 
and 3D. Our model, which does not make any explicit 
thermodynamical assumptions nor requires any fitting 
parameters, was capable of reproducing with great ac- 
curacy the results of recent experiments of the ATRAP 
collaboration Q. 
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